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ABSTRACT 

Employing data collected during the first 25 months' observations by the Fermi -LAT, we describe and 
subsequently seek to model the very high energy (> 300 MeV) emission from the central few parsecs of our 
Galaxy. We analyse, in particular, the morphological, spectral and temporal characteristics of the central source, 
1FGL J 1745. 6-2900. Remarkably, the data show a clear, statistically significant signal at energies above 10 
GeV, where the Fermi -LAT has an excellent angular resolution comparable to the angular resolution of HESS 
at TeV energies. This not only reduces dramatically the contamination both from the diffuse background and 
the nearby gamma-ray sources, but also makes meaningful the joint analysis of the Fermi and HESS data. 
Our analysis does not show statistically significant variability of 1FGL J 1745. 6-2900. Using the combination 
of Fermi data on 1FGL J1745.6-2900 and HESS data on the coincident, TeV source HESS J1745-290, we 
show that the spectrum of the central 7-ray source is inflected with a relatively steep spectral region matching 
between the flatter spectrum found at both low and high energies. We seek to model the gamma-ray production 
in the inner 10 pc of the Galaxy and examine, in particular, cosmic ray (CR) proton propagation scenarios 
that reproduce the observed spectrum of the central source. We show that a model that instantiates a transition 
from diffusive propagation of the CR protons at low energy to almost rectilinear propagation at high energies 
(given a reasonable energy-dependence of the assumed diffusion coefficient) can well explain the spectral 
phenomenology. In general, however, we find considerable degeneracy between different parameter choices 
which will only be broken with the addition of morphological information that 7-ray telescopes cannot deliver 
given current angular resolution limits. We argue that a future analysis done in combination with higher- 
resolution radio continuum data holds out the promise of breaking this degeneracy. 
Subject headings: Galaxy: center — synchrotron radiation: cosmic rays — molecular clouds: general 
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1. INTRODUCTION 

Over the past decade-and-a-half since the discovery by 
EGRET of a very high energy (VHE) gamma-ray source near 
the Galactic center (GC), there has been intense speculation 
as to what mechanism(s) are producing the observed emis- 
sion. The subsequent discovery of TeV gamma-ray emission 
from the Sgr A* region by the ground-based gamma-ray in- 
struments, in particula r by the HESS array of atmospheric 
Cherenkov telescopes (I Aharonian et al .1120041) . has generated 
further theoretical activity. Of general interest and import - 
given the GC constitutes the nearest example of a galactic nu- 
cleus - is the question concerning the sites and mechanism(s) 
by which particles are accelerated to TeV energies and beyond 
in the dynamical center of our Galaxy. 

Despite the fact that the GC TeV gamma-ray source is 
a point-like object for HESS, the 0.07° PSF of the instru- 
ment and the extremely crowded and complex nature of 
the region (as ev idenced by the complex radio morphology 
(lLaw et al.1 120081) ^ do not allow the unambiguous identifi- 
cation of the source(s) of gamma-ray emission. With the 
latest data, however, it is possible to place the center-of- 
gravity of th e TeV point sourc e within the central ~ 6" of 
the Galaxy (lAcero et al.1 120101) . leaving only a handful of 
possible sour ces. These include the central black hole it - 
self, Sgr A* (lAharonian & Neronovll2005al iLiu et al.ll2006T) : 
a plerion d iscovered within t he central few arcseconds of 
the Galaxy (TWang et al.ll2006t iHinton & Aharonian! [2007b . a 
putative "black hole plerion" pr oduced by the wind from 
Sgr A* dAtovan & Dermerl 120041) . and the diffuse < 10 pc 



region surrounding Sgr A* (lAharonian & Neronovl [2005b; 
Balla ntyne et alJl2007Ll20"Toh . 

Given the above background, we consider here the further 
insights now possible in light of the Fermi -LAT observations 
of the GC region. In particular, since the PSF of Fermi above 
~ 10 GeV is similar to that of HESS, it is possible to explore 
a quite broad energy interval of relativistic particles localized 
in this region. 

First 1 1 months of Fermi observations of the GC were pre- 
sented by J.Cohen-Tanugi on behalf of Fermi LAT collabo- 
ration during the 2009 Fermi Symposium. In this work, the 
authors argued that the Fermi source 1FGL J1745. 6-2900 and 
the HESS source J1745-290 are spatially coincident. Also, 
they derived an energy spectrum of the Fermi source till 100 
GeV and concluded that to match the HESS spectrum either a 
high-energy break or a cut-off is required. 

In this work we analyze 25 months of Fermi LAT data. In 
addition to the central GeV source and other reported sources, 
our analysis reveals four new sources of GeV gamma-rays lo- 
cated in this region. With the spectral information from both 
Fermi and HESS in hand, we model below the production of 
gamma-rays from the inner GC due to hadronic interactions 
of protons accelerated within the central black hole and dif- 
fusing into the surrounding interstellar medium. 

In section [2] we describe the reduction and analysis of the 
Fermi data. We present the details of our model in section [3] 
In section@]we discuss the implications of the obtained results 
and summarize the main conclusions in section[5] 
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2. DATA ANALYSIS AND RESULTS. 
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TABLE 1 

Coordinates and TS of the new sources discovered during 

THE ANALYSIS. 



RA 


Dec 


TS 


(J2000.0) 


(J2000.0) 




264.906 


-28.555 


331 


266.210 


-30.360 


424 


270.060 


-30.091 


189 


270.697 


-30.626 


192 



The large area telescope (LAT) on board the Fermi satellite 
is a pair-conversion gamma-ray detector operating between 
20 MeV and 300 GeV. The LAT has a wide field of view of 
^2.4 sr at 1 GeV, and observes the entire sky every two orbits 
(~3 hr for a Fermi orbit at an altitu de of ~565 km; full details 
of the instrumentation are given in lAtwood et al.l (|2009)). The 
data used for our analysis are based on the first 25 months of 
observations (August 4, 2008 - August 18, 2010). 

The data analysis was performed using the LAT Sci- 
ence Tools package w ith the P6_V3 post -launch instru- 
ment response function (Ran do & et al.ll2009h . The standard 
event selection for source analysis, resulting in the strongest 
background-rejection power (diffuse event class) was applied. 
In addition, photons coming from zenith angles larger than 
105° were rejected to reduce the background from gamma 
rays produced in the atmosphere of the Earth. The analysis 
was further restricted to the energy range above 100 MeV, 
because below this energy effective area becomes very small 
and the residual uncertainty in the instrumental response is 
significant. 

In order to take into account the broad point spread func- 
tion (PSF) at low (~ 100 MeV) energies, we constructed a 
sequence of test statistic (TS) images of the 10° x 10° re- 
gion around the Sgr A*. In producing TS images, we used the 
gttsmap tool with a tolerance parameter of ftol = 10~ 5 
and a bin size in each map of 0.1°. Finally, after subtract- 
ing the 19 known sources from the one year Fermi cata- 
logue (1FGL) which happen to be within the selected region, 
we found four new sources, which are listed in the Table [T] 
in the residual images. One of these sources (indicated as 
J1744. 8-3021) - shown in magenta in Figure Q]- lie within 
the 1.5° x 3° area around the GC. This source coincides spa- 
tially with known HESS source HESS J1745-303 and EGRET 
source 3EGJ 1744-301 1. 

In order to construct a light-curve for 1FGL J1745. 6-2900, 
we used a spectral method by selecting data in 300 MeV- 
100 GeV energy range and fitting all known sources, selected 
as above, with a single power law model. Afterwards, we split 
the whole time interval into 25 equal time bins and fit source 
spectra by fixing their slopes to the best-fit value obtained over 
the entire time period, leaving the source normalization as a 
free parameter. The normalization of the Galactic and extra- 
galactic background was also left as a free parameter. The 
resulting light-curve is shown in Figure [2] and is relatively 
stable and does not show any statistically significant varia- 
tion. The averaged flux is equal to (324.9 ± 7.05) x 10~ 9 
counts cm~ 2 s -1 , with a reduced \ 2 = LI f° r 24 degrees of 
freedom. 

Spectral fitting was performed within 100 MeV-300 GeV 
energy range with the gtlike tool. The spectrum in 
100 MeV-300 GeV energy range can be fitted by a power 
law with a slope of V — 2.212 ± 0.005 and a flux normal- 
ization of F = (1.39 ± 0.02) x 10~ 8 cm" 2 s" 1 MeV" 1 at 
100 MeV. We also attempted to split the spectrum into two 



different energy bands, and found that the fitted slope is equal 
to T = 2.196 ± 0.001 in 300 MeV-5 GeV energy range, and 
T = 2.681 ± 0.003 in 5 - 100 GeV energy range. The errors 
given above are statistical errors and represent the lcr devi- 
ation. Thus the slope of the Fermi spectrum above several 
GeV is significantly steeper than the spectrum reported by 
the HESS collaborati on at TeV energies (-FjjeSS ~ » 
dAharo nian et al. 2009). Note that at low energies, Fermi has 
a very broad PSF, rapidly moving from 4° at 100 MeV to 2° 
at 300 MeV. Thus, taking into account the possible source 
confusion in the region, one should treat the first point in the 
spectrum (100-300 MeV) with caution. 

3. MODELING 

As proposed in lAharonian & Neronovl d200 5b). a signifi- 
cant fraction of the protons accelerated near the black hole 
may enter the surrounding gaseous environment and initiate 
VHE gamma-ray emission through neutral pion production 
and subsequent decay. The efficiency of the process, and the 
energy spectrum of resulting gamma-rays depends not only 
on the protons' injection rate and the ambient gas density, 
but also on t he speed of proton transport into the surround- 
ing medium (lAharonian & Neronovl 1200551 iBallantv ne et al.l 
120071 120101) . To explain the gamma-ray spectrum r eported 
by the HESS collaboration, Aha ronian & Neronovl (l2005bl) 
assumed that relativistic protons with a power-law spectrum 
possessing a spectral index of V ~ 2 are injected into the 
dense gaseous environment surrounding the central black 
hole. The diffusion coefficient, D, was assumed to have 
a power-law dependence on energy of the form D(E) = 
10 28 (i?/lGeV)' 9 K cm 2 s _1 . For the cosmic ray diffusion 
in the Galactic disk k ~ 1 and f3 ~ 1, but of course the 
diffusion coefficient in the GC could be quite different. In 
lAharonian & Neronovl (|2005b) the paramete r (3 was assumed 
to be in the r ange of 0.5-1. In the model of Ballantvn e et al.l 
(2007, 2010) the propagation is treated using the ray-tracing 
technique. They found that in order to reproduce in their 
model the reported energy distribution of TeV gamma-rays, 
the spectrum of protons should be hard with a spectral index 
~ 0.75. Such an exceptionally hard injection spectrum of pro- 
tons implies a very strong energy dependence of the character 
of propagation of protons which, within the formalism of dif- 
fusion, would require a diffusion coefficient with (3 ~ 1.5. 

Given that the VHE emission detected by HESS and Fermi 
can be localized to within the central several arcminutes then, 
for a GC distance of d ~ 8 kpc, the linear size of the pro- 
duction region of VHE gamma-rays can be as large as 10 pc. 
Continuum X-ray and radio observations of the central 10 pc 
region of the Gala xy show a complex environment with many 
unique structures (lLaw et al.ll2008t ICrocker et al.ll2010b . For 
simplicity in this paper we assume that the supermassive black 
hole in the center of our Galaxy is surrounded by a shell of 
a dense matter with a density, normalized to nn = 1000 
cm~ 3 at 1 pc radius, with either constant or 1/r 2 radial de- 
pendence. The inner and outer radii of this shell are pa- 
rameters in our model. Another relevant parameter is the 
time evolution of the proton injection. Although one can 
treat it as a quasi-stationary process, in fact the proton in- 
jection can be dominated by one or several flares that oc- 
curred in the past in Sgr A*. In this context, one should men- 
tion the morphological interpretation of the diffuse gamma- 
ray emission observed by HESS from the central 200 pc re- 
gion of GC, which relates the positive detections of gamma- 
rays from giant molecular clouds in GC to a putative "proton" 
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FIG. 1 .— TS maps of the central part (1.5° X 3.5°) of the Galaxy center as seen by Fermi in 300 MeV - 3GeV, 3GeV-30 GeV and 3OGeV-30O GeV energy 
ranges (left to right). Positions of new sources are marked with magenta circles. Green ellipses correspond to the positions of the sources from the 1 year Fermi 
catalogue. Note that linear colour scheme has different maximum value in all cases varying from 5500 in the less energetic left picture to 140 in the most energetic 
right one. Source significance can be approximately estimated as a square root of TS. 

In the standard diffusion approximation the propaga- 

450 tion of particles is describe d by the diffusion equation 

dGinzburg & Svrovatskiil 1964b which, in the spherically sym- 
metric case, reduces to the form: 
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FIG. 2.— Lightcurve of the 1FGL J1745.6-2900 in 300 MeV 
energy range. The average flux is shown with a dashed line. 
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where n(r, t, E) is the space density of relativistic particles 
with energy E, at instant t being a distance r from the source; 
6(e) = —dE/dt is the continuous energy loss rate; Q(E,t) 
is the injection rate; and D(E) is the energy-dependent dif- 
fusion coefficient. We have assumed here, for simplicity, that 
D is independent of r and has a power-law dependence on 
energy as stated above. The solution of equation ((TJ can be 
written as (SyrovatsldHHH: 



flar e that occurred in Sgr A * in the past, 10,000 years ago or 
so dAharonian et al.l f2006). The detection of reflected X-ray 
emission from the Sgr B2 cloud is another, more direct piece 
of evidence about th e short flaring activity of Sgr A* a few 
hundred years ago dSunvaev et alJl!993l;lKov ama et ali i 19961 
l2008tlRevnivtsev et al.ll2004tlTerrier et al.ll2010t) .~ 



n(E,r,t)= [ P(E,r,t ~ x)Q{E,x)dx, 
Jo 

where the propagator, P(E, r, t), is defined as: 



P(E,r,t) 



[4^A(S,i)]3/2 



exp 



A\(E,t) 



(2) 



(3) 
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and 



X(E,t) = - \ 

JEJ, 



E g (t) &0) 



(4) 



In equation E g is the energy that a cooled particle has at 
time t, if its initial energy was E. 

Formally, the diffusion equation does not contain infor- 
mation on how fast a particle may propagate. Since Eq.(fl]i 
does not prevent an artificial "superluminal motion" (w = 
2D(E)/r > c), we follow the p henomenological approach 
proposed bv lAloisio et al.l (120091) who introduced a propaga- 
tor, P(E, r, t) in the form: 



P(E,r,t) 



1 



a(E,£) 



47r(rf) 3 (1 -£ 2 ) 2 Ki[a(E,Z)\ 



exp 



(5) 

where 6{x) is the Heaviside step function, £(t) = r/ct, K\{x) 
is modified Bessel function of the second kind, and a(E, t) is 
defined as: 

a ^ = 20jy (6) 

In the low-energy regime (i.e., E <C E c ), the propagator 
given by Eq.© reproduces the standard treatment of diffu- 
sion, whilst in the high-energy regime (E > E c ) it describes 
particles that move in a rectilinear fashion. Here E c is the 
energy at which Eq. (HJ allows diffusion with the speed of 
light 



E r 



( cR 



\2Dr 



GeV. 



(7) 



Due to the energy dependence of the diffusion coefficient, 
proton propagation will be quite different at low as compared 
to high energies. We have explored how this plays out in 
the GC environment. The result presented in Figure [3] shows 
the change of the radial distribution of protons as a func- 
tion of energy. We have determined the proton distribution 
after 300 years of continuous injection into the interstellar 
medium of density n# = 10 3 cm~ 3 within a region of ra- 
dius R — 3pc. The initial spectrum of protons was assumed 
to have a power-law distribution with an exponential cutoff, 
Q(E) oc E- 2 exp(-£/100 TeV). 

It can be clearly seen from Figure [3] that, whilst at 10 TeV 
the particles pass through the region in an almost rectilinear 
fashion, at lower energies, protons propagate diffusively. An 
explanation for the lack of low-energy protons at high radii 
can be found in the fact that these particles have too low an 
escape velocity to travel that far in a given time. 

The spectrum of the protons integrated over the gamma-ray 
production region (see Figure [3]l is shown in Figure [4] Pho- 
tons produced by the interaction of relativistic protons with 
such an energy distribution fit both Fermi and HESS data well. 
At low (GeV) energies, the diffusion radius is smaller than the 
region so that protons are accumulated within the region and, 
given the almost energy-independent pp cross-section, mirror 
the spectrum of the injected protons. On the other hand, at 
TeV energies protons begin to propagate in a rectilinear mode 
and will have again the form of the injected spectrum, albeit at 
a lower flux level. Protons with an intermediate energy have a 
much steeper, diffusion-processed spectrum representing the 
transition between the two extremes. The spectral shape of the 
highest energy gamma-rays is not affected by the propagation 
effects. Therefore in order to match the spectrum at highest 
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FIG. 3. — The fluxes of protons as a function of radius at different energies 
shown as marked. For ease of comparison each energy was multiplied by 
a factor of 900 (for 10 GeV), 8000 (100 GeV), 5 X 10 6 (ITeV) and 7 X 
10 s (10 TeV). The flux at each energy have been multiplied by B? , so that 
rectilinear propagation corresponds to a horizontal line. 
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FIG. 4. — The energy distribution of protons averaged over 3pc of gamma- 
ray production region, as reconstructed from Fermi and HESS data. 
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FIG. 5. — Spectral energy distribution of gamma-rays expected from a re- 
gion filled with relativistic and non-relativistic protons within different as- 
sumptions concerning the injection, diffusion and the region geometry (see 
text for a discussion of parameters for each specific model). The data points 
have been derived from the Fermi and HESS data 



energies reported by HESS (Aha ronian et al.l 120091) . we as- 
sumed an exponential cut-off in the proton spectrum and fix 
its position at 100 TeV. 
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Below, we fit parameters which represent the particle in- 
jection spectrum, the propagation of the injected protons, and 
the geometry of the interstellar medium. It is instructive to 
systematically examine the influence of these model param- 
eters on the resulting spectrum. To do this we begin with 
10 4 years of injection of relativistic protons with a spectrum 
of the form Q(E) oc £r 1!) exp(-.E/100 TeV) into a 3 pc 
radius region filled by an interstellar gas of constant den- 
sity, riH — 10 3 cm . The injection rate was taken equal 
to Qo = 3.9 x 10 39 erg/sec. The diffusion parameters were 
chosen as f3 = 0.9, k — 0.015 in order to reproduce the com- 
bined Fermi and HESS spectrum. The photon spectrum re- 
sulting from this parameter set (model A) is shown in Figure 
|5] The other curves in Figure [5] illustrate the effect of change 
of a single parameter, while all others are fixed to the values 
used in model A. 

Figure[5]corresponds to the case of a source active for only 
300 years (model B). This change does not affect high energy 
particles, traveling rectilinearly, because their escape time is 
t esc = R/c ~ 10 years, much shorter then the injection time. 
These particles fully fill the region and their density is the 
same as in the case of model A. The diffusion time at low 
energies (E < E c ~1 TeV in this case), however, is much 
longer (viz. t diff = R 2 /2D ~ 9000 years at lGeV). Thus, 
300 years will be not enough for low-energy particles to travel 
to the outer regions of the shell, and the total spectrum is ex- 
pected to be harder with respect to model A. Thus, in this 
parameter set, whilst the radiation does not differ from our 
heuristic case at high energies, at lower energies there are nec- 
essarily fewer gamma-rays. 

Models C and D in Figure [5] show what occurs when the 
diffusion parameters k and (3 are changed. If one increases 
the diffusion coefficient by a factor of 10, then by the same 
factor, the diffusion time of the low-energy particles is de- 
creased, leading to a corresponding reduction in the intensity 
of the gamma-ray emission. If one changes the energy de- 
pendence of the diffusion by decreasing (3 by a factor of two, 
then the transition of the particles propagation from diffusion 
to rectilinear propagation occurs at much higher energies (E c 
has 1 / f3 dependance, see equation (0, and reaches PeV en- 
ergy in this case). Thus the emission increases at all energies 
and the spectral form changes due to a larger influence of high 
energy particles. 

Models E and F in Figure [5] show the change of the photon 
spectrum due to the change of the geometry of the region. The 
corresponding curve in figure corresponds to a shell geometry 
with an outer radius of R = 10 pc (as compared to the 3 pc ra- 
dius considered in the other models). Using this particular ge- 
ometry causes the overall normalization of the resulting spec- 
trum to increase (a factor of 0.25 was applied to this spectrum 
for easier comparison in the production of Figure [5]) due to 
the increase in size of the gamma-ray production region. Ad- 
ditionally, the 'bump' at lower energies becomes wider, since 
the energy at which particles transition from the diffuse to the 
rectilinear propagation regimes increases as R 1 ^, and thus 
the larger the radius, the larger the number of 'intermediate', 
GeV, energy particles that can be accumulated. 

The effect of the shell volume is clearly seen in the case 
when the value of the inner radius of the disk is changed. 
Model F in Figure [5] is for a shell with inner radius R m i n — 
2pc and outer radius R max = 3pc. In this case the overall 
normalization of the resulting spectrum decreases (in Figure 
[5] a factor of 3.3 has been applied to this spectrum in order 
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FIG. 6 — Brightness profile of the GC in lOOMeV - lGeV (solid lines) and 
lGeV - lOGeV (dashed lines) energy ranges. See Table|2]for the description 
of the models. 25" corresponds to lpc. 



TABLE 2 

Parameters of different models used for Figure[6] The 
medium density at r=lpc was assumed to be equal to 
n h = 10 3 cm -3 . qo and qf correspond to the injection rates 
of the constant source and a flare correspondingly. 



model 


QH 


K 


P 


Qo 

10 39 erg/s 


pc 


Q/ 

10 39 erg/s 


const 


const 


0.02 


0.95 
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1/r 2 


1/r 2 


0.002 


0.75 


0.1 








flare 


const 


0.7 


0.5 
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8000 


shell 


const 


0.02 


0.95 
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to more easily aid the comparison to other models), but the 
low-energy bump is still wider than in the case of a shell with 
no hole (i.e. an oblate spheroid). This is due to the fact that 
removing the inner part of the shell mostly diminishes the soft 
photon emission, effectively increasing the relative number of 
more energetic ones. 

Finally note that the sharp drop of gamma-ray spectra below 
1 GeV is the result of the kinematics of pion production at p- 
p interactions, and therefore does not depend on the model 
parameters. 

The radial distribution of photons is also highly dependent 
on the model parameters. Figure [6] shows the brightness pro- 
file of the inner 3 pc after 10 4 years of the constant injec- 
tion at energies between 100 MeV-1 GeV (solid lines) and 
1 GeV-10 GeV (dashed lines). For all models the constant 
source was active for 10 4 years. The initial spectrum of pro- 
tons is assumed to be a power-law with an exponential cutoff, 
Q(E) oc iT 2 exp(-£yi00 TeV). The medium density at 
R=lpc was assumed to be equal to nji = 10 3 cm -3 At higher 
energies particles pass through the region almost rectilinearly 
and the region will appear as a point-like source of gamma- 
rays to Fermi . As labeled in Figure|6l different curves corre- 
spond to different model parameter sets in terms of the ambi- 
ent matter distribution and injection. For all models in Figure 
|6j parameters, given in Table|2] were chosen so that the result- 
ing integrated emission accurately reproduces the Fermi and 
HESS data, and the resulted profiles were normalized to the 
maximum value to aid comparison. 

Geometrically, the solid line represents the case of constant 
density and exhibits a broader profile at higher energies (line). 
If instead, one models the region with a density falling off 
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proportional to a tih oc 1/r 2 profile, the resulting profiles 
are thinner and are represented by the red solid and dashed 
lines, which almost coincide on the figure. The profile is more 
centrally peaked in the latter case, since the matter is more 
concentrated in the center and so the photon flux will originate 
mostly from this region. 

The green solid and dashed lines show the profile created if, 
in addition to a constant source, there was also a flare which 
occurred 300 years ago (the injection rate in the Table|2]cor- 
responds to a flare of length 10 years). In order to match ob- 
servational data in this case, a larger diffusion coefficient had 
to be assumed, which inevitability leads to a larger diffusion 
radius and - correspondingly - to a wider profile. Finally in 
Figure [6] we show the brightness profiles corresponding to a 
shell geometry. In this case the profile has a maximum at a 
radius of the inner shell. 

4. DISCUSSION 

The spectral properties of the very high energy emission 
from the GC differs considerably from that at lower energies: 
whilst the GC is known to be variable at X-rays and near-IR 
wavelengths, no variability has been detected either by HESS, 
or by Fermi. This seeming duality has a natural explanation if 
the low energy emission is generated very close to the central 
black hole, while the gamma-ray emission originates from a 
much larger region and is emitted during the diffusion of the 
relativistic protons through the interstellar medium surround- 
ing the central black hole. In such a case, the very high energy 
emission would only reflect (with a delay) major flares orig- 
inating from the central source. As remarked above, one in- 
terpretation of the distribution of the diffuse, TeV 7-ray emis- 
sion relative to the molecular clouds in the central ~ 200 pc 
of our Galaxy, is that a central CR pro ton source flared about 
10 4 years ago dAharonia n et alJl2006l) . In the previous sec- 
tion we showed that our model is able to reconstruct data 
if we assume a constant injection of relativistic protons for 
1.0 x 10 4 years (see Figure|5]l. This time is higher than the 
diffusion time for the used set of parameters and thus the ob- 
tained photon energy spectrum is effectively steady state. 

The observations of reflected X-radiation from the cloud 
Sgr B2 located at a distance of 100 pc from Sgr A* suggests 
that a few hundred years ago there was an i ncrease of X- 
ray luminosity of Sgr A* (|Sunvaev et al. 1993; Kovam aet al.l 
[T996ll200irRevnivtsev et al.l2004tlTerrier et alj2010t) . In our 
modeling, we checked that we are able to explain the data as 
a result of a constant injection of protons over three hundreds 
years. We found that if one assumes an injection rate of pro- 
tons of 6 x 10 39 erg/s, and takes /? = 0.95, n = 0.01 (the 
resulting radial distribution of protons for this model is shown 
in Figure[3]l, then the resulting emission is in good agreement 
with the observations. Our model is able to self-consistently 
explain different spectral indices at GeV and TeV energies by 
the different effective escape velocities of the protons. While 
high energy protons, producing TeV photons, escape quasi- 
rectilinearly without spectral deformation, as, indeed, do the 
particles fully trapped at the lowest energies, the particles with 
intermediate energies are affected by diffusion, but not fully 
trapped and their spectrum becomes much steeper providing 
the transition between two extreme cases. 

Recent monitoring of the Sgr B2 cloud with X-ray in- 
struments shows fl ux variability on time scales of 10 years 
dTerrier et al.ll2010l) . This variability can be naturally inter- 
preted as a measure of the flare duration. Although the X-ray 
luminosity and proton acceleration in Sgr A* should not nec- 
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FIG. 7. — Combined Fermi (green points) and HESS (blue points) explained 
by superposition (black solid line) of a proton flare of 10 years duration 
happened 300 years ago (magenta dashed line) and a constant source that 
switched on 10 4 years ago (red dotted line). See Section|4]for model param- 
eters. 

essarily correlate, it is interesting to explore also the scenario 
when we deal with a flare of proton acceleration and injection 
into the interstellar medium on timescales of years. 

In Figure|7]we compare the spectra of gamma-ray emission 
resulting from realization of three different scenarios: (i) a 
proton flare of 10 years duration that occurred 300 years ago, 
(ii) a constant source that switched on 10 4 years ago, and (iii) 
a proton flare on top of the constant source, namely the super- 
position of (i) and (ii). To fit the data, we took the size of the 
gamma-ray emission region to be R = 8 pc, parameters of 
the diffusion coefficient j3 = 0.65, k = 1, and initial proton 
spectrum of Q(E) oc E~ 2 cxp(-i?/100 TeV), with the pro- 
ton injection rate of 1.9 x 10 39 erg/s for the constant source 
and 1.9 x 10 42 erg/s for the flare. 

For this parameter set, the 300 year old flare cannot have 
a strong impact on the observed TeV spectrum at this point 
in time, since most of the high energy protons from the flare 
have already escaped. On the other hand, the emission at GeV 
energies is produced by protons from the flare which are still 
trapped by diffusion in the gamma-ray production region. To 
explain the TeV data we need much slower diffusion or a fresh 
injection of protons, for example contributed from a very re- 
cent flare, or by the quasi-steady component of protons. Ac- 
tually the form of TeV emission doesn't depend on the age of 
the source if it exceeds t esc — R/c ~ 30years. The case of 
superposition (solid line) of the flare (dashed line) and persis- 
tent (dotted line) components of protons is shown in Figure[7] 
For the chosen parameters, the GeV energy range of gamma- 
rays is dominated by the flare component of protons, while 
the TeV gamma-rays are contributed mainly by protons from 
the persistent component. 

Thus, we are able to reproduce the observed broad-band 
spectrum of gamma-rays in different ways. In all cases the 
required injection rate is well below the Eddington limit of 

the 4 x 10 6 A/ black hole in GC, L Edd = 47rGMmp ~ 

q ' C( j T 

6 x 10 44 erg/s. The total energy required in relativistic protons 
currently trapped in the gamma-ray production region varies 
from 10 49 to 10 51 erg for different models. This energy can 
be injected in very different ways: in reality there has proba- 
bly been a series of flares with different energetic signatures 
occurring throughout the life-time of the central source. 

The observed spectral and temporal properties of the GC at 
various wavebands are not enough to constraint all the param- 
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eters in our model. Additional information can be extracted, 
in principle, from the gamma-ray morphology of the inner 
arcminute region: above, we showed that different set of pa- 
rameters describing the observed spectral properties result in 
very different radial profiles (e.g., Figure [6). However, with 
the angular resolutions of the current space- and ground-based 
detectors, we cannot distinguish between the different radial 
profiles. Fortunately, such information can be recovered by 
observations of synchrotron emission of secondary electrons 
from decays of charged pions, accompanying the production 
of gamma-rays from decays of neutral muons. Since through 
this channel the electrons and gamma-rays are produced with 
similar energy distributions, we can connect directly the fre- 
quency of synchrotron photons of secondary electrons with 
the energy of the "genetically" connected gamma-rays: 

^ 100 (h^g)(iSv) 2mHz ' (8) 

where the magnetic field is norma lized to the probable mini- 
mum value expected in the region ( Crocker et al.l l2010). 

Thus, in the first approximation, the morphology of the syn- 
chrotron radiation of "hadronic" origin should be similar to 
the morphology of GeV gamma-rays. While at sub-GHz fre- 
quencies, GC radio photons are attenuated by free-free ab- 
sorp tion in dense Hll re gions between the GC and the Earth 
(see lCrocker et al.ll2010l) . the synchrotron emission at ~ GHz 
frequencies arrives without significant attenuation. The emis- 
sivity of synchrotron radiation of secondary el ectrons in the 
inner f ew pc of the Galaxy has been studied by lCrocker et al.1 
(120071) , based on a model where the interactions of pro- 
tons, diffusing away from an assume d central source, suppl y 
the observed, point-like TeV signal (B allantvne et al.ll2.Q07l) . 
They then compared the predicted synchrotron emission to 
the GHz radio frequency spectrum, and found that essentially 
all the diffuse non-thermal GHz radio emission from the cen- 
tral parsecs of the Galaxy could be explained as due to emis- 
sion of secondary electrons (and positrons). Therefore we an- 
ticipate that the new Fermi data combined with available ra- 
dio measurements, should allow us to constrain significantly 
the parameter space of models positing that the GeV and TeV 
gamma-ray emission of the GC is due to hadronic interactions 
in the central few parsecs of GC. Analysis of the morphology 
of radio emission holds out particular promise here. The re- 
sults of such an analysis are beyond the scope of this paper 
and will be presented elsewhere. 

It should be noted that cosmic ray electrons produce pho- 
tons not only in the radio domain; their bremsstrahlung emis- 
sion can be an important source of high energy gamma- 
rays. Comparing the synchrotron and bremsstrahlung cooling 
times: 

W-4xl0-(^) 2 (^) \ (9) 

^™- 1 - 2xl ° 12 (loc^)" ls ' (10) 

we conclude that bremsstrahlung losses dominates at energies 



below: 

This suggests that bremsstrahlung will dominate at energies 
less than ~ 30 GeV for fiducial nu and B values. 

The relative importance of electron bremsstrahlung in pro- 
ducing the observed gamma-ray emission, is characterized by 
the ratio of cooling times of electrons and protons associated 
with the bremsstrahlung and neutral pion production, respec- 
tively: (ftf/q* ~ (3t pp /t br )f ~ 4/, wh ere / = n e /n p i s 
the electron to proton energy density ratio dAharonianll2004l) . 
Hence, if the ratio of protons to electrons is 3> 1, then ir°- 
decay gamma-rays dominate over bremsstrahlung. The con- 
tribution of electron bremsstrahlung to TeV gamma-rays is 
quite limited because of the severe energy losses of very high 
energy electrons due to the synchrotron and IC losses. 

Finally, we note that the Fermi data pres ented here can not 
be exp laine d by IC models prop osed by lAtoyan & Dermerl 
(2004) and iHinton & Aharonianl (12007b . While these mod- 
els are in a good agreement with HESS data, they predict that 
the energy flux in the GeV part of the spectrum should be 
smaller than in TeV range, apparently contrary to the Fermi 
observations. 

5. SUMMARY 

We have analyzed 25 months of Fermi data on the GC re- 
gion. The Fermi LAT source 1FGL J1745. 6-2900 lies within 
the error box of HESS source J1745-290. We found that, 
while below 5 GeV, the spectrum of 1FGL J1745.6-2900 has 
a photon spectral index similar to the HESS source, the spec- 
trum at higher energies is better described by a steeper spec- 
tral index. We have formulated a model which produces 
a photon spectrum that can naturally explain the observed 
broad-band gamma-ray emission. This model considers the 
hadronic interactions of relativistic protons which, having dif- 
fused away from a central source, presumably the central 
black hole, fill the inner few parsecs of our Galaxy. We have 
explored the parameter space of our model, in terms of the 
geometry, characteristics of the diffusion coefficient, and in- 
jection rate history. 

We have shown that the available spectral information can 
be well described with different sets of model parameters and 
that additional information is required to distinguish model 
scenarios. Such information could be obtained from the spa- 
tial distribution of the observed gamma-ray emission; how- 
ever, the required arc-second resolution cannot be reached by 
gamma-ray telescopes. Luckily, synchrotron emission from 
the secondary electrons and positrons expected in our model 
may be detected by radio telescopes which possess an angu- 
lar resolution high enough for the purposes of distinguishing 
between model parameters. 
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cussions on the transition from diffusion to rectlinear regime. 
The authors wish to acknowledge the SFI/HEA Irish Centre 
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Ireland (SFI) under its Research Frontiers Programme. 
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